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ABSTRACT 


Miewescritatery torces generated by the interaction of 
eeralie@: gravity waves with a rigid horizontal circular 
cylinder located near a simulated ocean floor were deter- 
Maned experimentally. A circular cylinder of six inch dia- 
meter was subjected to gravity waves of wave lengths varying 
from 2 to 18 feet in a two-dimensional wave channel. Both 
the horizontal and vertical components of the hydrodynamic 
force were measured. 

Experimental results are plotted in dimensionless 


Boe the horizontal and vertical dimensionless force co- 


Memeterents are plotted as a function of the relative water 


L. 2 .hais fo ae Deak (er Gere ie 3 PL eae te, mney gl ye ae 
PeesscraeilvecewWwave LCP LI, alu Cie PoLadtive wave 


@epei, 
Moment. lt was tound that the horizontal force coefficient 
mantied quite tinearly with the wave height, and the maximum 
weve Was approximately equal in either direction. fThe 
vertical force component, however, showed a much more com- 
plex variation. The upward and downward magnitudes were 
approximately equal for small amplitude waves, but as the 
wave height increased, the variations became quite non- 


linear and the magnitudes in the two directions differed 


Considerably. 
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NOMENCLATURE 


Deep Lon 


Cimenirarecy Linder 


added-mass coefficient, dimensionless 


depmeneon submergence of cylinder 


maximum horizontal force per unit length 


FL /[ega* (H/2a)], hoetzonitaletoree coefficient. 


dimensionless 


Maximum positive vertical force per unit length, 


(upward) 


F./ [ega? (H/2a)], POsteive Vertical foree, Ccoctti- 


Creve dimenstonless 


Leinui iecitavesvemtlcal f£oree per Unit length, 


(downward) 


Fy/Tega* (H/2a)], negative vertical force coeffi- 


cient, dimensionless 


qccclewarion due to cravity 


ne eh Go. 


the incident wave 


water depth 


wave number= 27/L 


wave length of the incident wave 


emetn Ot 
perted, Of 
Crime 

component 
component 
component 


component 


eile ey leimnicen 


the incident wave 


Sreclimamveclocttyetimetne X-dtrect ion 
eee liemaece Heme rolein Ehe x-direction 
SMe lOcity In they direct von 


Cieepmlnidwaccele pation, 1m the y-direction 





Symbol PeSecripet Lon 


x pee COO natid ce 
y Space coordinate 


MeloOcilty potential 


0 densormeyeor tne fiudad 
O frequency of the incident wave= 217/T 
wT Mia hee id Viscosity 





d/a 
F/[oga* (H/2a) ] 
F/ [ega* (H/2a) 
H/2a 
ny a 


gna] L 


relative 


DIMENSIONLESS NOTATION 


depth of submergence 


hoplzonltal wave force coet ficient 


vert aca | 
Gelacilye 
relative 


relative 


wave force coefficient 
wave height 
Wate pedep tn 


wave length 
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I. INTRODUCTION 


MimeEcGcwimn, canst nere has been an accelerated interest 
in the exploration of the ocean bottom with special atten- 
men devoted to those depths encompassed by the Continental 
emelt. this interest has been motivated in great part by 
the desire to exploit the mineral assets of the ocean bot- 
tom. Additionally, investigation of ocean bottom pipelines 
memecopecitically applied to sewage outfall lines, cil pipe- 
Iimes and large aqueducts are under consideration. 

As a result, there has been considerable interest in 
the hydrodynamic interaction of gravity waves with large 
Seemarine structures resting on the ocean bottom. The 
design of large submarine offshore structures necessarily 
heduines One to have an understanding of the forces induced 
697 such environmental phenomena as surface waves. To this 
pom the theoretical work has yet to be a conclusive tool 
Midssisting the design engineer to predict all expected 
wnarerons to whieh the anchored structure will be subjected 
as it lies on the ocean oor. Experimental observations 
Gietne forces induced by ocean waves are difficult to ac- 
Guerre, and 1t 1s theretore necessary to simulate as pre- 
eisely as possible the conditions existing at the ocean 
bottom. 

Extensive experimental work has been carried out since 
Stokes, who in one of the earliest papers on the subject 


[1], showed that, by experimentation with a pendulum-type 


Jha 





Preclmtacmexyhessiom £LOr the force on a cylinder oscilla- 
merci ail intinate fluid consists of two terms, one in- 
volving Meee at on Smee cylinder. and the other 
mvolving velocity. 

Wave forces exerted on piles have been intensively in- 
vestigated during the past two decades. It can be shown 
Mumeewinen the relative size of the body, i.e., the cylinder 
meeameter, 15 small compared to wave length, a condition 
Pomciets generally Satisfied in the case of piles, the 
familiar Morison equation (aletsea Valid approach. this 
eiieatron combines both a drag term and an inertia term, 
and the values of the Spemimmenteabl cOcrtriclents of drag 
Biamimertia are left to the ingenuity of the designer. At 
Pose "muemomemmenoncal coctficients arc dependent 
upon the geometry of the objects under consideration and 
the amplitude of the fluid motion, and it would appear that 
420 years of research has shed little light on understanding 
Emmis santerdependence. Moreover, the proximity of the free 
Surface as well as the impermcable bottom surface affect 
tncs & coefficients. ie Vic; Ommenese wait leuluies. thie 
Morison equation is not relied upon in the present approach 
to the problem, and the results are represented directly in 
dimensionless form rather than through the drag and inertia 
eOert 1 CLenet . 

Keulegan and Carpenter [3] used periodic standing 
waves in their analysis of a cylinder subjected to a field 


Peccecmiubatorny motton. they were able to separate the total 
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ao wotraicemintomdrapy and inertial components. They correla- 
ted values of inertia and drag coefficients under conditions 
Sivarying Gurrent intensity and cylinder sizes with a peri- 
Oa parameter, U_T/D, a dimensionless group combining the 
maximum fluid velocity U> Ene mwave perivoa I, and the cylin- 
Gemeaiameteér. This parameter is a form of the relative 
immiid displacement since UL? woepEOnporte1onal te the ampli - 
fuee or the fluid motion. 

Sarpkaya and Garrison [4] studied both theoretically 
ome xperimentally unidirectional flow past a circular cyl- 
meer with Constant acceleration. They also showed that 
the total resistance as well as the drag and inertia com- 
Mements Of force may be represented as a function of rela- 
@evenr+iluid displacement. In the present analysis of the 
experimental data, this same general kind of dependence on 
mind aisplacement 1S supposed, and, accordingly, the force 
coefficients are correlated with the ratio of the wave 
Menent to cylinder diameter. 

JohmsoOlmep so feadseanvestigated the horizontal forces 
resulting from the interaction of gravity waves with a 
wimewlan cy lanaer placed relatively near the bottom. since 
only long waves were studied, the resulting forces were 
supposed to be dependent on the water depth and wave height 
emt.) this study shed little light, however, on the basic 
problen.. 

Heubuchmandsohank. 1m a recent study [6], have con- 


—~ 


ducted wave channel testing of vertical and horizontal loads 


aise 


Om submerged three-dimensional structures due to interaction 
with gravity waves. Several models were studied using a 
two-dimensional wave tank, and they concluded that the re- 
meets Of these experimental studies were sufficient to be 
epplicd to preliminary designs. Unfortunately, few actual 
experimental points were presented and therefore the degree 
of correlation is somewhat suspect. 

Saimersomnee: gl. (7), investigated the wave/structure 
interaction problem as applicable to larger objects where 
miemori1son equation 1s no longer valid and a diffraction 
Bimeery Must be employed. In this situation, the incident 
wave is mericcred DpOhmeneGolnteh~meuthe large structure, 
and the assumption that the object does not affect the in- 
mmcmuevavre To mo loncer valid. The diftfracticn 
meeounits tor the relative size of the body and, neglecting 
Paecous SeecOUsmmccUMesu tne cCxlstence Of a velocity poten- 
Mitte wince Laree structures, such as the recently deployed 
Simeeaso Bridge and Iron Persian Gulf oil storage tank, are 
woarened, the analysis must include consideration for rela- 
mervensize and £ree surface effect, and this basic consid- 
ekation 2S inherent in the diffraction ue 

Lemappearsrenat there will continue to be heavy re- 
Pianeceron experimental model results in the design process 
until more reliable theoretical predictions of wave induced 
forces are assured. To this time, large structures have 


been placed on the ocean floor only after extensive model 


muvGctesuiave been Carried out. 
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Micmpresenussumay was undertaken to obtain an insight 
Mirco tne interaction Of gravity waves with large submarine 
meanuetcires Located at Cuailecmerle ocean £1loor.) The study 
meemudes direct measurement of horizontal and vertical 
merces On the submerged cylinder. The results are pre- 
femecd in dimensionless form as a function of the relative 
feeeer depth, the relative wave length, and the relative 
wave height. 

iiemincerest in attempting this investigation was 
generated while participating in numerous discussions on 


Ziemeroulem wumder the guidance of Dr. C. J. Garrison. 


Ws 


Meee wieOret TCAL CONSTDERATIONS 


Mec cieiatke nepresentation of the problem under con- 
Mmi@eretion 1S Shown in Figure (1). A rigid circular cylin- 
MmemeOot radius a submerged with center at depth d in water 
memaepth h 16 subjected to a train of regular waves of wave 
Pemeeiele aia herent H progressing in the positive x-direc- 
mow lie primary interest lies in the horizontal and ver- 
Piealecomponents of force induced by the waves on the 
Simeecular cvlinder. 

Riacmeonplele mmeoretie¢al Solution to the present prob- 
ies ditfticult, even assuming the flow to be unseparated; 
and, therefore, in order to proceed with an experimental 
meres tausaupncecessary tC anply dimensional anal- 
ese to the problem. It is known @ prtort that the maximuin 
(monizental or vertical) force per unit length of cylinder 


Pmcependent upon the following variables: 


eo” P ymax = £(h,d,a,L,H,o,g,u) - (1) 


where h denotes the water depth, d the depth of submergence, 
pecne cylinder radius, L the wave length, H the wave height, 
op the fluid density, g the acceleration of gravity, and u 
Gicmelund viscosity. The period of the incident wave is not 
required in the above relationship because a second rela- 
tionship already exists between the wave period and the three 


vortantes fo H; and L. [his is true not only in the case of 
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— 


x 
H 2 
d 
q cylinder 


X \ | channel floor 


hIotincwlens D6Finition Sketen . 


sy, 


linear waves, but in the case of large amplitude non-linear 
waves. 

A dimensional analysis of the physical quantities in- 
fmedted in Equation (1) yields the following dimensionless 


Variables: 


F F 
ae eee e/a. d/a, H/2a, —= 


pga* pga* L OVgha? ] 
Mees last term on the right hand side of Equation (2) rep- 
Mmeeemes the ratio of the Froude number and the Reynolds 
ieee reand, accordingly, is an indicator of the ratio of 
mmeous f0rces to gravity induced forces. If this number 
Meelarge, «corresponding to a large scale flow and fluids 
G@eeecmall wiscositv, it may be stsnected that the gravita- 
meomal £Orces may dominate the flow, and the viscous ef- 


HeerssmMay be Considered unimportant. As a result, this 


parameter inay be disregarded and Equation (2) may be writ- 


ten as: 
Y scaieEss P max 21a 
——— or oie ——— ya. d/ a, H/ 2a GS) 
pea piceae lb 


A second reason for the disregarding of this parameter 
may also be cited. In unidirectional motion starting from 
rest and moving with constant acceleration, it has been 
eho bye oarpkaya and Garrison [4] that the total force 
Sxenecd, oma Gircular ner can be expressed as a func- 
miOneGtetheenolatrive displacement of the fluid only. Al- 


though the present problem involves oscillatory as opposed 
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Pomniloirectional Motion, it is nevertheless suspected that 
the same general kind of dependence on relative displace- 
ime would resuit. In the present problem, the amplitude 
Oeeene relative fluid motion is proportional to the wave 
Mement to cylinder diameter ratio and, therefore, H/2a 
Mommesents an equivalent parameter. In the experimental 
meectam, an attempt is made to correlate the two force co- 
efficients on the left hand side of Equation (3) with the 
dimensionless ratios on the right hand side. 

Meenoucnetne complete theoretical solution to the pro- 
posed problem is very complex, it is possible to obtain an 
approximate solution if certain simplifying assumptions are 
mooted. § Ihe velocity potential associated with the regular 
MopressiVe WaVe dione 1s given Cy: 


cosh k(yth) 
cosh kh 


=) 
i 
SR 


cos (kx-ot) (4) 


amd strom this expression, the velocity and acceleration of 
wicmiiuad oarticles within a wave may be written, respec- 
tvely; as: 


gk> cosh k(yeh) 


u =" cosh kh sin (kx-ot) (3) 


O 
gk sinh k(yth) 
O 


Ve == cos (kx-ot) (6) 
oki 

u = ——A cosh Kien) cos (kx-ot) (7) 

: eke sinh k(yth) 

Vv = a ee ote (8) 





O cosh kh 
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In wave force problems involving objects whose dimen- 
sions are small in comparison to the wave length of the in- 
SmdeTit Wave, it 1S Common to assume that the object does 
Meceattfect the incident wave. It is then assumed that the 
Geyect 15 subjected to a uniform oscillatory flow having 
acceleration equal to that induced by the wave at the cen- 
Memmomecne Object af the object were not there. Following 
mess Line Of reasoning, the horizontal force per unit length 
meemcylinder may be written as the product of the virtual 


mreseand ftluid acceleration as 


i = oma? (1+C_) Uma, (9) 


if it can be assumed that flow separation does not occur. 
BUMetrbuting LOY UW From Byuation (7) and evaluating at 
Seoeecaves the dimensionless horizontal force coefficient 


as 


IK 


271 
cosh 2™ (h-d) | 
spar = n(1+C,) 272 | ——EE | i :. (10) 


cosh ath 


Since the dimensionless force coefficient is linear in the 
wave height, this expression may be written more conveniently 


as 


cosh oa (h/a - d/a) 


F 
x 27a 
= 1+C ee 11 
pga? (H/2a) e | y | cosh At2 h/a _ 


The parameter oe jeclurring Ine Bowatrom (i) denotes tie 


added mass coefficient and depends on the shape of the object 
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meguestion as well as the proximity of the bottom and free 
euortace. Values for the configuration in question which 
corresponds to a circular cylinder between a rigid boundary 
ana free stirface are not available to this writer's know- 
Mert linc attempc im this thesis, therefore, is to deter- 
meme by trial and error the value of Cs VG ee yc Sar etaic 


mest fit to the experimental data. 


Ze 








Peewee iteNtAlL INVESTIGATION 


‘A. APPARATUS 

The wave channel as depicted in Figure (2) was de- 
eeoned with a rectangular cross section two feet deep with 
metath of 15 inches. The channel was fabricated in sec- 
tions. Each section was composed of two vertical walls 
micmome Norizontal floor member. The vertical walls were 
pteeroned by use of three 4 foot long 2 x 4's attached to 
eeamen xX 8 foot x 3 foot sheets of plywood at the ends 
and middle of its longest dimension. The channel floor 
ilember was composed of two 15 inch x 8 foot long plywood 
Sheets separated by a 1 inch x 4 inch strip around the 
poem ine bottom and side members were mated together 
Wiechietong 5/8 inch threaded rods and formed the basic 
Giammeimecross section. Six of these 8 foot sections were 
mrcimpolted together to form the basic channel. sSquare- 
Mess of the channel cross section was maintained by use 
of spacers placed between the stiffening members near the 
mhoor. Every butt joint was given a heavy application of 
Silicone sealant to prevent leakage. 

At a location 30 feet from the wave generator end of 
the channel, a 3 foot wide plexiglass observation window 
extending the full height of the channel was installed in 
both walls of the channel as shown in Figure (3). The 
test cylinder was then mounted in the center of this win- 


dow ‘ 
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Figure 3. Wave Channel Observation Window. 


24 





Dissipation of the generated waves was provided by the 
peOomottUcwlnemiocated im the last section of the channel. 
The beach as shown in Figure (4) consisted of three layers 
Grmperttondted Stainless steel sheeting. Layer separation 
acme ructural framing were provided by 2 x 4 runners with 
Premmmbareest cross SE€ction dimension perpendicular to the 
Sheeting. The slope of the channel beach was maintained 
Ateron | . 

The wave generating component of the channel as shown 
in the photograph of Figure (5) had two basic elements, the 
piston-type plate and the driving motor. The plate was an 
aluminum sheet 1/8 inch thick with adjustable edges on the 
Vemmmedi sides @lving 1t the capability of adjusting to 
GOMmsorm tC ce siaes Oo: the chanmel. The pilatc was mounted 
on a framework attached to 4 slider bearings, and the 
Sede mepeiarings rode on two 120 inch brass rods attached 
LOmercmenanmele walls eeA drivange rod was attached to the 
Diudwemstructuralhemembers and connected to a flywheel at- 
Licieamugmene OuUCpUtL Shatt Ofva Vari-Drive motor. The 
(mitemnocmids attached towthe flywheel by means of a bolt 
CidmelOtetm such a Manner that the Ree aneniel ty could pe 
varied between zero and 6 inches. In this way, waves of 
various amplitudes could be generated. The variable speed 
transmission was driven by a 2 hp inductance motor and pro- 
Vidieimmapccdmnuaice cw tenaime «trom ZZ to 180 rpm. giving 
Veco Mien rromels r£eet to 1.5 feet at a water 


depth of 18 inches. 
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Pepeexigle@ss Ci1rcular cylinder of 6 inch diameter with 
an end clearance of approximately 1/16 inch was used in the 
best reported herein. As depicted in Figure (6), the cylin- 
ger was hung by two fine wires at the appropriate distance 
off the wave channel floor. Weights were added inside the 
meander to pretension the wires and to prevent slack in 
eieewires during testing. 

Piicmuppcmmends Of the wires were attached to small 
alumanum Cantilever beams equipped with strain gages on 
Deen the Compression and tension side of the beams. These 
webtieal tOree Strain pages were then connected into a 
Wheatstone bridge and read-out on one channel ae Wwe = 
channel Sanborn Model 296 amplifier/recorder. The cylinder 
Mame cti aimed in wie horizontal 2irectzion by similar fine 
wires connected into a harness and around a 4 inch diameter 
ball bearing pulley mounted on the wave channel floor as 
shown in Figure (7). These wires were then attached to a 
Dai Ot Cantilever beams equipped with strain gages, and 
the Strain gases were connected into A Wheatstone bridge 
nlG@mtoad-OuEsOnmeLne second Channel of the two-channel am- 
plairer/recorder . 

The four beams, two for the vertical force and two for 
the horizontal force, were attached to a plywood mounting 
board positioned over the channel as shown in Figure (8). 
Hie wOcamanwebtomall exactly alike, and the section of the 
beams where the strain gages were mounted was reduced to a 


emeossmesection of 1/4 x 1/4 inches from a nominal beam cross 
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Figure 7. Horizontal Wire Harness and Pulley Assembly. 
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Figure 8. Beam Mounting Board. 
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Seetron ot 1/2 x 3/8 inches. The machining provided a 
SiC@ummEmItcdclOneInto the area of stress concentration. 
oe -aeA EZ oot Zo spay gages were iplred toluene top 
PAGmbOtteniesutiaces Of the beams. Each beam was 9 inches 
long and had mounting holes for two moment arm positions, 
SHemaeno Inches sand one at S anches from the strain gage 
MOGtelCiew this variable load attachment point gave each 
beam a sensitive force measurement capability under a 
itaecmmancec Orevalues, Ihe beams were designed to with- 
Stand a six-pound load at the largest moment arm. 

The wave height gage, as shown in Figure (9), con- 
sisted of a small diameter aluminum Fone c ine as the 
idstemenaie to which a plexiglass mounting table was at- 
facwede ss Une elements of the Wheatstone bridge were 
mounted on this table, and two 30-gage copper wires spaced 
1/2 inch apart were extended from the table and attached 
to an insulator mounted at the end of the aluminum rod. 
An 18-inch separation between the levels of the table and 
insulator provides a large measurement range capability. 
The other end of the rod was mated to a traverse block 
stand which allowed the gage to travel + 6 inches, and 
provided a convenient, yet accurate, means for calibration 
Dy eadjmsctinowene submergence of the wires. The complete 
bitt was posacaoned over the channel and centered at ap- 
proximately 5 feet in front of the cylinder. 

The wave height probe was of the parallel wire re- 


ml-tanee type. This type of gage opérates on the principle 


Sz 








Figure 9. Wave Height Gage. 





Pade tie Conductance between the two parallel submerged 
wires varies approximately proportionately to the length 
Siewine Submerged. The wave gage was connected in parallel 
momoue ley of a Wheatstone bridge as depicted in Figure 
wry. 10 compensate for the inability of the carrier am- 
Memerer tO Offset the resistance of the immersed parallel 
foeees, an e€xternal variable J0K ohm potentiometer was con- 
Memecd aAcrOSS af Opposite leg of the bridge. This resistor 
was used for coarse balancing of the bridge. The Wheatstone 
Gimdse Was connected via signal and excitation leads to a 
emetuamle "Garrier preamplifier and read out on a one-channel 
fee Sanborn ROMA Ionececrder: 

Mimincnone limntiary Operation while attempting to clear 
EitemproObplems £Lrom the system, it was Wound that installation 
of an impermeable wall between cylinder and channel floor 
Was a necessity. The experiment was designed to simulate 
treeplacement Of the circular cylinder on the bottom, but 
1t was necessary to maintain a small but finite gap between 
Sviinder —anad the bed of the channel to insure that the 
pair of wires connected to the cantilever beams provided 
the sole vertical support. However, it was found that the 
Memeical torces as well as, to a lesser extent, the iO Wie 
zontal forces were strongly dependent on the gap width. 
Ehowevisualigation studies showed that a high velocity 
NeemOceianucd aamethne small gap under the cylinder. Thus, 
Mimowdewstomelinimate the flow under the cylinder and yet 
not influence the force measurement system, a flexible 


plastic sheet was installed across the gap between the 
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meemncer and tloor. ihe sheet was joined to the cylinder 
by pressing the plastic into a machined slot in the cylin- 
feoevall> sand ethen anchoring it with a length of "0O"'-ring. 
Mie other end of the plastic was laid under a metal sheet 
meeanches x 14 1/2 inches and 1/16 inch thick and placed 
ae the channel bettom as shown in Figure (11). The gap 
Demvicen thew cy lander and tleor was adjusted to approxi- 
mately 1/8 inch. This gap was found to be large enough 
Poi sSuremunctecne cylinder meyer touched the bottom when 
Slloyectedmio etme Largest wave. 

lise order PomeaciOlNocmtaromunecenaener large intluence of 
Poss sialsmcapmundenm the Cylinder, Figure (12) was prepared 
showing both the maximum horizontal and vertical force as 
a function of wave amplituue. [It iS apparent from this 
fesne Cidteamemovibwon thes plastic barrier had the effect 
of reducing the maximum horizontal force a slight amount. 
With the barrier in place, the upward force was generally 
quite large while the downward force was rather small. 
However, upon removal of the Dopitcmmrarather high veloc- 
acy tlowetimeuchmiiessapmwas Observed by injection of dye 
mear thewwave channel floor, and this een region of high 
velocity apparently resulted in low pressures on the under- 
Side Of@eicwe: Minder cCausmc ,ehe downward force to be in- 
Greased secomemacuaoiye Ine Lasher sextreme effect of the 
gap is cléarly demonstrated by comparison of the various 


eurves indicated on Figure (12). 
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Inasmuch as the intent of the present study was to 
determine wave forces on a cylinder lying on the bottom 
mierein tie gap would not be present, the plastic barrier 
was maintained in place throughout the test program for 


eieeeeases where the cylinder was placed near the bottom. 


Be feol PROCEDURE 

iii formers ecep in the test procedure was to fill the 
Peewee nannelswith water to obtain the desired water depth 
namermimdersdiameter ratio, h/a, Experiments were con- 
ducted at the three depths, 12, 15, and 18 inches. After 
Memccired Water level was attained; a careful inspection 
Waieeedtried OU tO insure that all equipments were properly 
“prepared HOT Oper delon- 

Pinemviel LL Cai i rempOstdulelM serews were adjusted to 
einem ocsSinedwelevattvom of the cylinder off the chan- 
Hem DOteem. The horizontal Witocmmere Dretensioned to 
Pacey Cieoreie king and the attendant errors im force measure- 
ments. 

A typical run may be described as follows. With water 
adjusted to the proper depth, a given wave length, or equiv- 
alently, wave period was established by setting the speed 
Cael orto rive tilt. subsequent to a sufficient warm-up 
period, the amplifier/recorders shown in Figure (13) were 
balanced and readied for calibration. Calibration of the 
three recording channels, horizontal force, vertical force, 


and wave height was carried out at the attenuation expected 
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Pomel USC diuming the run. Aber calibration, the wave 
height was set by adjusting PiemceGcontrrcity Of the wave 
generator to the appropriate amplitude. The wave generator 
was then turned on, and after a steady state wave was es- 
tablished in the channel, the recorders were turned on and 
data was taken. The wave generator was then switched off, 
pGCmeemicxuewave height was set on the flywheel eccentric 
VigefemMdiitainming the same speed setting. Wave heights at 
the fixed wave length were varied within the capability of 
Ciemeyowcn tO produce répresentative readings usually ranging 
mPommaporoximately 1/4 inch to 5 inches. At the conclusion 
Of Wech Series of Poicmteeeivene wave length, @ check cali- 
UEotomron thie three channels was quickly carried out. Then, 
Gumtewmwave Length was Selected and Set Of the Vari-Drive, 
mdmeremcyele of balancing, Calibration, wave height selec- 
tlopmmmand wave production was repeated. 

Calibration was performed on the load cells as follows. 
The recording chart was set at 1 /sec. Speed ,eand Loads in 
Qeamemidans temsuwere placed on the saddles located beneath 
and hung from EicmioneZzontal Sensing wires as shown in Fig- 
temo cCaltbratiom loads were placed initially on the 
after beam, then on the forward beam. The vertical load 
cell was calibrated by use of a platform of known weight 
Ccllecmead OMmememeylinder as shown in Figure (15). Weights 
Wer eCmenen Glaeedeon this platform. 

Wave height gage calibration was performed by adiusting 
the vertical position of the gage through the use of 4 slid 


bearing and screw. The procedure was carried out as follows. 
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Horizontal Force Calibration Operation. 


Figure 14. 
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The gage was raised 1/2 inch, simulating a decrease in water 
level, then lowered 1/2 inch below the null, simulating an 
Mieaease Iiewater levels Ihen, the gage was raised +1 inch 
emombewered -ll imch, etc. This procedure closely simulated 
the actual passage of a wave profile with crests and troughs. 
Mie cenerated wave crest passing through the test location 


is shown in Figure (11). 
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ieee Cte NAL RESULTS 


As indicated in Equation (3) in the THEORETICAL CON- 
mebbRAlIONS, the wave force coefficients can be represented 
aemerunctions of relative water depth, relative depth of sub- 
Mmemmecnce, relative wave length, and relative wave height pa- 
rameters. Accordingly, a test program was devised to measure 
miemwave £OLrce coefficients for practical ranges of the in- 
Mependent parameters. Four series of tests were conducted 
at different water depths and depths of submergence. Three 
memenese Corresponding to depth to zylinder radius ratios 
Poo. 0, and 6.0, were Carried out with the cylinder 


iawmene HOCLOM, 1,e€., at d/ja= 3.0, 4.0, and 5.0, respec- 


mee iy ) A fourth series of tests was conducted with the 


Cylinder raised one radii off the bottom corresponding to 


(yee 4,0 at a water depth ratio of 65.0. 


Xe Pye uC Aim WwAVE. FORCE TRACES 

Lemont mnr ero mdihtrclile tO completely characterize 
tme force Variations with a single number, Cree) Ser tala 8 
GiimemGead=Olts mor the Various series of tests are pre- 
Seco eiilee douress(to-28). Figures (16-19) show the hori- 
Zoomed and Vertical force Variations with the cylinder on 
CicmIOEEONTCOrhesmolaine to a relative depth of submergence 
of d/a = 5.0 and a relative water depth ratio of h/a= 6.0. 
The figures are arranged according to increasing values of 


the wave length parameter ranging from 27a/L= 0.09 to 0.97. 
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On each figure, data corresponding to four values of rela- 
tive wave height are presented in ascending order. 

Figure (16) representing the largest wave length and 
having a depth to wave length ratio of h/L= 0.08 corre- 
sponds to waves which are somewhat outside the sinusoidal 
range. These long waves are characterized by long, smooth 
Mmeeuens and Sharp crests. Accordingly, the horizontal 
force traces are characterized by positive and negative 
peaks which result from the maximum acceleration occurring 
Momemunechestae Of the crest. ~fhese two peaks are followed 
by a relatively Pit MOheLOneNGitathe torce trace Corres - 
monmaine tO the long, smooth trough portion of the wave. As 
wane Seen trom the figure, this kind of variation appears 
memo redOmiliale Even as tne wave neighi 1s increuseu for 
this long wave case. 

ieee verttredietorce, however, shows a much more incon- 
mistent variation. If the system natural frequency is dis- 
regarded, definite upward peaks are observed at points on 
the trace corresponding to the wave crests. These peaks 
ere apparently related to the reduced pressure over the 
Mpper portion of the cylinder which is associated with the 
velocity squared term in the Bernoulli equation since the 
maximum velocity occurs at the crests. In sinusoidal waves, 
an equal and opposite velocity occurs at the troughs which 
euils Naturally give rise to a second upward peak in the 
force trace. However, in shallow water waves, the veloc- 
iemeassOchateauwithrha trough is much reduced and, there- 


fore peaks corresponding to the crests only are observed. 
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Pitimemyscorrcsponds to a shorter wave length having 
imme bativeswave length value of 27a/L= 0.16. It is apparent 
Prommcic ttelre that both the horizontal and vertical forces 
vary sinusoidally for the smallest amplitude wave correspond- 
mieeto H/2a- 0.5. As the wave height is increased, it is ap- 
Mearent that the horizontal force continues to vary sinusoid- 
aumieye tor the full ramge of amplitudes. The vertical force, 
however, shows a sinusoidal variation for the smallest am- 
plitude waves and becomes increasingly less sinusoidal as 
maecewave height is increased. Moreover, two separate peaks 
Mepear to develop on the vertical force trace as wave height 
ancreases, one associated with the maximum velocity occur- 
Mmneeate cme crest and the other associated with the trough. 
Hecomaximum velocity at both tWe crest and the trough re- 
Zabeoerm Gecuiced pressures om Che upper portion of the cylin- 
femmcid, Consequently, corresponds to upward vertical forces. 

Figures (18) and (19) correspond to still smaller wave 
Memeths of 2nva/L= 0.24 and 0.77, respectively. It may be 
noted that again the horizontal force shows a sinusoidal 
Variation throughout the range of wave amplitudes while the 
vertical force shows the sinusoidal variation at small am- 
plitudes becoming non-sinusoidal as the amplitude is in- 
creased. Figure (19), however, which corresponds to the 
shortest wave lengths shows sinusoidal variations through- 
out the range of wave amplitudes for the horizontal as well 


aemevertucal. force. 
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Figure 17. Wave Force Trace, nh/a= 6.0, d/a= 5.0, 2v7a/L= 0.16. 
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Horizontal and vertical force coefficients Corres ponud-— 
ing to relative water depth of h/a= 5.0 with the cylinder 
@m@ethe bottom, i1.e., d/a= 4.0 are shown in Figures (20-22). 
imtesce three figures represent three different values of 
melative wave length extending from 21a/L= 0.10 to 0.67. 
iimesseneral trends tend to be similar to those corresponding 
to the larger depth of h/a as shown in Figures (16-19) and 
previously discussed. Data corresponding to the smallest 
value of 2rTa/L as presented in Figure (20) shows trends 
typical of nonlinear shallow water waves. The horizontal 
force shows the null period associated with the long, smooth 
Meoughs, while the vertical force is characterized by the 
upward force peak corresponding to the wave crests. As the 
Pre men eels acGreascad, Che Berizontal force acain dis- 
ferays Sinusoidal variation as indicated in Figures (21) and 
(22). For the shorter wave lengths, the vertical force shows 
@esimusoidal yariation at small amplitudes with superimposed 
ipgward peaks OCCUrring near the maximum velocity location at 
Mee crests and troughs. Figure (22) corresponding to the 
Shortest wave length shows a sinusoidal variation for both 
mnree Components throughout the range of wave heights. 

ine tr OCee Vd tattom wecords for the series of tests 
eomeduavcdmauan,d=—4.0 are presented in Figures (23-25). 

Three representative values of relative wave length, 21a/L= 
en wana s45eare presented, Figure (23) shows varia- 
tions similar to those of equal values of 2na/L at the two 


Fate ‘nl ae 


larger depths already discussed. Figures (24) and (25) 
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Bo esiMibrar tO previous runs, the primary difference being 
mie increased magnitude occurring at the smaller depths. 
iicmainmatinserNes Of runs was Carried out at a relative 
water depth ratio of h/a= 6.0 with the cylinder suspended 
eieeradii Off the bottom, i.e., d/a= 4.0. Typical force 
meteeseat threc ditferent values of relative wave length, 
may b= 0.10, 0.16, and 0.6, are presented in Figures (26), 
mae, and (28), respectively. The horizontal force trace 
as shown on the figures reflects the same general trends as 
Mmremonreviousty described series. The vertical component 
memrorce also shows trends similar to those previously dis- 
cussed with phe: @yinder on the bottom, but the peaks as- 
meoGrated with maximum velocity in the wave occurring at the 
rests and trouphs are seem tu be in the opposite direction. 
iireers, the vetocity of the fluid through the one radii gap 
Between the cylinder and the floor is Miporenil yy COs 1de¢ r— 
ably larger than the velocity through the space above the 
cylinder and beneath the free surface and, therefore, a 
mMcking-dowh d@clion aS Opposed to a lifting action results. 
Pip uevMOuse Tums, liemay be moted that Figure (28), cor- 
responding to the largest value of 2na/L, shows a sinusoidal 
fertation tor both force components throughout the range of 


wave heights. 


Be PaVveEwoOrcreCOmerICIENTS 
In order to represent the results of the experiments in 


a simpler, albeit less descriptive form, the maximum values 
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of inewdamenmsitomless £orce Coefficients obtained from the 
force traces were plotted as a function of relative wave 
height, H/2a. Results corresponding to the relative water 
amereor n/a wo,0 with the cylinder on the bottom, i.e., 
dja= 5.0, are presented in Figures (29-32). Results cor- 
Pesponding to the other two relative water depths of h/a= 
Paeeand 4250) also with the cylinder on the bottom are pre- 
femecd 1 Fagures (33-36) and (37-40), respectively. These 
results correspond to values of 2na/L ranging from approxi- 
meaeety 27a/L= 0.08 to 0.97. 

As was apparent from the force traces, the horizontal 
force showed sinusoidal variations with equal magnitude in 
Mmemmedirections and, therefore, only one curve was required 
memecpresent ticeamplitude Of this Variation with relative 
wave height. The dimensionless force coefficients are de- 
fined as f=F/pga*, and three curves are noted on each fig- 
Breo,eene fOr the horizontal and one each for the upward and 
aogigard Vertical force. As is evident from the figures, 
BhemnOlrlzZontal foree Coefficient Shows a linear variation 
with the wave height in all cases. 

Howeviel aS etiewvertical torce showed, in general, larger 
upward ohen downward magnitudes, the two curves labeled . 
and f are presented showing the maximum upward and downward 
Mic omOtmene= Vehencal toree COctTlicient, respectively. 
Throughout the range of the data, the vertical force coeffi- 
cient for the cylinder located on the bottom showed a quite 


Sousiotentevardacion, the upward force being always larger 
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4 


than the downward force and having a generally concave up- 
ward trend. The downward force was always quite small show- 
mie trearly linearly increasing variation followed by a 
breaking olf at higher wave amplitudes. 

The variations of the force coefficients as observed 
Maeriegures (29-40) can be justified if one considers the 
components which make up the dynamic pressure acting on the 
Mimeace Of the cylinder. Two components of pressure are 
Menmcrally recognized, the one component being associated 
feos the Unsteady term in the Bernoulli equation, and the 
feeomd being associated with the velocity squared term oc- 
eur ring genet a. Pitesulcstead,y = temmarmvolyes the velocity 
Seetential to the first power and, therefore, is linear in 
H/Za. The second term involves the spatial dérivative o 
Eiemvelocity potential oe Wee ING iS eeeeWere lore, propor - 
mecmeMen to (i/2a)-. thus, linear variation of the horizon- 
@ileworce with H/2a apparently results from contributions 
PaOlmeteninsteady Component of pressure only. The velocity 
Squared term which would contribute a nonlinear component 
Omeeoree Call comtribute nothing to the horizontal force if 
tiem low Pattern 1S symmetrical about the y-axis. This, 
Gemectr oc, lo cexpeeted provided separation does not occur 
meamaGccOUuntomeonm tne linear Varlation in the horizontal 
force with the parameter H/2a. 

HileCIwmerEdledls force variation with H/2a, however, shows 
markedly different trends. Both the upward and downward 


i als 4 


components seem to vary linearly with relative wave 
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Figure 38. Wave Force Coefficients for h/a= 4.0, d/a= 5.0. 
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Bomesiall amplitude waves, and the upward force tends to 
bend upwards while the downward force tends to decrease at 
larger values of H/2a. These variations, however, are con- 
Seoeelt With the previous discussion regarding the pressure 
Simeenwe slintace of the cylinder. In the case of the vertical 
momce, the flow pattern 1s not symmetrical about a horizon- 
Beeaxis passing through the center of the cylinder, and, 
Bieterore, the contribution from the velocity squared term 
mimeeie Bernoulli equation does not cancel. As a crest or 
muon passes the cylinder, the wave induced velocity is a 
feinun Giving rise to a component of pressure (or force) 
proportional xo (Vee. thesvenruleal Troree coefficient 
Retnmatclon with H/2a, therefore, is composed of two terms, 
Mem Lopottional Uoll/Za and one proportional to (H/2a}*. 
Thus, its somewhat parabolic shape is not unexpected. 

item <inum values Of the horizontal and vertical force 
SeextiGicnts ror the tinal series of runs SOndule te deat ad 
water depth ratio of h/a= 6.0 with the cylinder suspended 
PemtieebOottom, 1.¢., dja= 4.0,"are presented in Figures 
fvemuenmough (44), As in the previously described tests, 
the horizontal force coefficients showed a linear variation 
with wave amplitude whereas the vertical force coefficients 
Showed an initial linear portion followed by a nonlinear por- 
Bion Onevenm ulna this Series Of runs, the high velocity 
in the gap beneath the cylinder eAlcedvaunet nonkinear,con- 
tribution to the force in the downward direction making the 
magnitude of fy larger than actond. Opposite tO Tome 


associated with the cylinder near the bottom. 
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Figure 42. Wave Force Coefficients for h/a= 6.0, d/a= 4.0. 
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c, HORTZO MEA FORCE COEFFICIENT 

The maximum values of the horizontal force coefficient 
ae Poeoniov meine noules(29-40)meare seen to vary linearly 
with the relative wave height, H/2a. In view of this lin- 
Prev eMavVOl wees POSSiIbDIe to represent these forces in 
a much more condensed form. That is, since the force co- 


efficients are linear with H/2a, the parameter, 


FL /[pga* (H/2a)], 


mmmecneactUally represents the slopes of these curves can 
Ipeepmesentcd ih a2 Single figure as a function of the rela- 
meee wave length, Zma/L. Such a plot is presented in Fig- 
Me w(4o) wherein this force coefficient 1s represented as 
Peeumetion Of Zina;L for tne four series of tests carried 
out at different water depths and depths of submergence. 
The upper three curves shown in Figure (45) represent 
three different water depths wherein the cylinder was lo- 
ated on the bottom. Nome kcCemoOumwater depth 1S clearly 
maicated, as the curves corresponding to decreasing depths 
lie above the one at a maximum water depth of h/a= 6.0. 
litte cenmeral trends of these three curves ao a leees Ite tee 
maximum values of force coefficients occurring near 21a/L= 
feo Wiel decreasing morce at shorter wave lengths. At 
longer wave lengths than those corresponding to the maximum 
value for the force coefficient, the curves tend. to become 
flat while the linear theory would indicate a decrease 


Similar to that associated with the curve corresponding tv 
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EXPERIMENTAL APPROXIMATE ANALYSIS 


BATS EQUATION (11) 
- N/a 
—--— 4.0 
eee OU) 


6.0 
eee 18 





4 

OZR CIS lS Clim O.2 O4 
Relative Wave Length, 27a/ | 

itso > Och Zontdl rorce Coctficients. 


83 


Oe Sls 0 





Aas Haven, await se possible that this behavior 
mmeaused by the tact that the wave is becoming a nonlinear 
shallow water wave. 

iiicenonuzZenmd et erGcemGoctiTc1ents corresponding to the 
Pourem series of tests carried out at a water depth ratio 
of h/a= 6.0 and with the cylinder suspended at one radii 
off the bottom is also shown in Figure (45) as the lowest 
curve. The general trends are seen to be somewhat similar 
to the cases where the cylinder was located on the bottom 
Beteept that the Curve tends to be somewhat flatter. 

DiecieaccenpenLOncolpare: the experimental results with 
mie approximate theory, Equation (11) was programmed on a 
G@ucital computer, and values for the horizontal force co- 
Pe nonesmni ere CCIeracCdugaG 2 sumetionm cf 2ra/lassuming 
various values for the added mass coefficient, cae These 
Powlesewerer then plotted and compared with the experimental 
results in order to determine which value of oe gave the 
Dewees INnesemeurves representing tlre “best fit’ were 
then plotted on Figure (45) and were compared with the ex- 
perimental plots. As might be expected, it was possible 
to obtain good agreement between the approximate theory 
and the experimental data for the largest depth ratios. In 
these cases, the proximity of the free surface, which was 
disregarded in the approximate analysis, was probably not 
Dart Lcularly amportant. However, OG thew shallowest depth 
where the free surface was more influential, it appeared 
to be impossible to find a value for C, which gave reason 


apllesaerecment with the plotted experimental results. 
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Ones poanteworeimeerest aim eigure (45) was to compare 
the two curves associated with a maximum depth of h/a= 6.0, 
mre One Case Corresponding, to the cylinder located on the 
bottom, d/a= 5.0, and the other with the cylinder elevated 
one radii off the bottom, d/a= 4.0. As can be seen from 
fae tigure, the latter curve 1s considerably below the for- 
Mer and commesponds @o an added mass coefficient of approxi- 
mately C= 1.2. The data corresponding to a similar 
eituation with the cylinder located on the bottom appears 
boragree relatively well with the curve based on ee = 2.8. 
thes marked increase in C., corresponding to the case where 
imencylinder was located on the bottom apparently results 
PPonmeene proximity of the bottom. 

A second interesting feature in Figure (45) was 
apparent discontinuity in the two curves associated with 
the two shallowest steps near 27a/L= 0.3. The cause of 
mere discontinuity 1S not certain, and it might be dis- 
mecca as Scatter in the experimental data. a iieieew is 
however, a second explanation which is equally probable. 
Geersepossipile to have a small amount of reflection at 
imc value ot 2na/lLeeand af 1S possible that this reflec- 
ted waves affected the wave height measurement. Depending 
on the phase of the reflected wave, it may either add or 
MbiitecteanaNecausSemal apparent increase Or decrease in the 
value of the wave force coefficient. The absence of this 
discontinuity im the curves corresponding to the greater 


depths is in agreement with this hypothesis since retlec- 


Bromeadecreases with relative depth. 
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V. CONCLUSIONS 


An experimental study of gravity wave interaction with 
aesubmerged horizontal cylinder has been conducted. From 
meme results Of this study, the following conclusions are 
drawn. 

Pe Cm newman Ge Gt Wave lliements tested in the experi- 
iment, the horizontal forces resulting from gravity wave 
interaction with the cylinder increased linearly with wave 
IMmemeoht . 

2. the vertical force was always much smaller than 
the horizontal force when the cylinder was located on the 
meottom. ihe downward part of this torce fluctuation was 
Pee AO (ee Cm ew tres ers Suowcd a nonlinear in- 
Grease with wave amplitude. At very smal] amplitudes 
faleree ene VelOeGity Squared Contribution to the pressure 
Was small, a sinusoidal variation of the vertical force 
Pesulted . 

Cee CU enemeacC GOI TCsponGing to the cylinder sus- 
pence One Tadiisort the bottom, the horizontal force was 
reduced and the downward force was greater than the upward 
force. 

Tweet OurGascececuuiere the depth Of submergence was large, 
the approximate analysis appeared to give valid results 
meovaded = the proper value of the added mass coefficient was 


iM OWwn.. 
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S. When the value of 27a/L becomes small, correspond- 
ing to long waves, the wave force trace shows features dis- 
Similar to those associated with the more sinusoidal waves. 
Bimparticular, a mull period in the horizontal wave trace 


Mecurs during the passage of the trough. 
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Experimental Wave and Force Data, h/a= 96,0, d/a— 5.0. 
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TABLE II - Experimental Wave and Force Data, h/a= 5.0, d/a= 4.0. 
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Experimental Wave and Force Data, h/a= 4.0, d/a= 3.0. 
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Experimental Wave and Force Data, h/a= 6.0, d/a= 4.0. 
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Wave Forces With and Without Barrier, h/a= 6.0, d/a= 5.0. 
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